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SUMMARY 


A modified form of Donnell* e equation for the equilibrium of 
thin cylindrical shells is derived which is equivalent to Donnell’s 
equation but has certain advantages in physical interpretation and 
In ease of solution* particularly in the case of shells having 
clamped edges. The solution of this modified equation by means of 
trigonometric series and its application to a number of problems 
concerned with the shear buckling stresses of cylindrical shells 
are discussed. The question of implicit boundary conditions alBO 
is considered. 


INTRODUCTION 


During a general theoretical investigation of the stability 
of curved sheet under load* a method of analysis was developed 
which appears to be simpler to apply than those in general use. 

The development of this method Is presented in two parts* of which 
reference 1 Is the first and the present paper the second. The 
specific problems solved by this new method are treated In detail in 
other papers. (See* for example* references 2 to 7.) 

In reference 1. the stability of a stressed cylindrical shell 
was analyzed in terms of Donnell's equation* a partial differential 
equation. for the radial displacement w* which takes into account 
the effects of the axial displacement u and the circumferential 
displacement v. Reference 1 shows the manner in which this equation 
can be used to obtain relatively easy solutions to a number of 
problems concerning the stability of cylindrical shells with simply 
supported edges. The results of the solution of this equation were 
shown to take on a simple form by the use of the parameter k 
(similar to the buckling— stress coefficients for flat plates) to 
represent the state of stress in the shell and the parameter Z 
to represent the dimensions of the shell* where Z is defined by 
the following equations: 
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For a cylinder of length L 

Z = l£ 1 1 - u 2 

rt ' 

and for a curved panel of width b 



where 

r radius of curvature 

t thickness of shell 

and 

p Poisson’s ratio for material 

The accuracy of Donnell’s equation was established by comparisons 
of the results found by its use with the results found by other 
methods and by experiment. 

In the simplest method that has been found for solving 
Donnell’s equation, the radial displacement w is represented 
by a trigonometric series expansion. This method can be used to 
great advantage for cylinders or curved panels with simply supported 
edges but leads to incorrect results when applied uncritically 
to cylinders or panels with clamped edges. 

In the present paper an equation is derived which is equivalent 
to Donnell’s equation but is adapted to solution for clamped as 
well as simply supported edges by means of trigonometric series. 

This modified equation retains the advantages of Donnell’s equation 
in ease of solution and simplicity of results. The solution of the 
modified equation by means of the Galerkin method is explained, and 
the results obtained by this approach in a number of problems 
concerned with the shear buckling stresses of cylindrical shellB 
are given in graphical form and discussed briefly. Boundary conditions 
implied by the method of solution of the modified equation are 
also discussed. 
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SYMBOLS 


a 

h 

a iJ> 

k-B 


P 

r 

t 

u 

v 

w 

X 

y 


length of curved panel (longer dimension) 
width of curved panel (shorter dimension) 


deflection coefficients in trigonometric series 

for cylinder; 


shear— stress coefficient 
,2 


rtb 


for curved panel 





dlrect-axiad -stress coefficient 




■■■ for curved panel t 

y 


{oxtL 2 

Vlht 2 


for cylinder; 


D 


/ +r d. 

circumferential-stress coefficient ( — 7 ■ * • • for cylinder 

\ V^Dir^ 

2 ~ for curved panel) 


D 


Dir 

lateral pressure 


/ 


radius of curvature of cylindrical shell 

thickness of cylindrical shell 

displacement in axial (x— ) direction of point on shell 
median surface 

displacement in circumferential (y— ) direction of point 
on shell median surface 

displacement in radial direction of point on shell 
median surface; positive outward 

axial coordinate 


circumferential coordinate 


i, 3, aCl 

n, P, <LJ 
D 


integers 

plate flexural stiffness per unit length 


/ Et 3 \ 

\ 12(1 - n 2 )J 



k 
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E 

F 

L 

Q» Qi> 
Rs 

Ex 

Z 


\ 

P 

ox 





Young’s modulus of elasticity 


Airy’s stress function for median— surface stresses produced 

/'£2 f 

( — s, stress in axial direction; 

Vdar 

— - j stress in circumferential direction; - r 

&c 2 . Sac 3y 


by buckle deformation 


shear 


stress 
length of cylinder 
mathematical operators 

shear-stress ratio; ratio of shear stress present to 
critical shear stress when no other stress is acting 

axial-compressive— stress ratio; ratio of direct axial 
stress present to critical compressive stress when 
no other stress is acting 

curvature parameter ( — i/l - p 2 for cylinder; 

\rt ' 


bf 

rt 


/l — p, 2 for curved panel or long curved strip) 


half wave length of buckles; measured circumferentially 
in cylinders and axially in long curved strips 

Poisson’s ratio for material 

applied axial stress , positive for compression 
applied circumferential stress, positive for compression 
applied shear stress 
critical shear stress 

operator + _ + a— . 


r 
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w 




V'^ 


inverse operator defined by equations 
= V^(V-^f) = f 


THEORY 

Derivation of Modified Equation 


The equation of equilibrium for a flat plate may be written 


IfA* 


+ t(cr 


x _ 


b^v 

bxP- 


+ 2r 


b^v 
2tx by 




= O 


( 1 ) 


where p is lateral pressure. (This equation is equivalent to 
equation ( 197 ) of reference 8.) 


For a cylindrically curved plate having a radius of curvature r, 
the following pair of simultaneous equations of equilibrium may be 
written (as a generalization of equations (11) and (10) of reference 9) : 


w\r 



+ 2t 


b^w 
5x 5y 



+ t d% = 

r d x 2 


0 


( 2 ) 


9^f — S. = 0 
r 8x 2 


( 3 ) 


where F is Airy r s Btress function for the median— surface stresses 
produced by the buckle deformation (reference 10) . Equation (2) 

t d^F 

r c* 2 ' 


differs from equation (1) only in the addition of the term 


which expresses the_effect of the curvature. Equation (3) shows 
that,, unlike flat plates 3 cylindrical shells experience stretching 
of the median surface when originally straight lines in the surface 


V 
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are bent slightly. Elimination of F between equations (2) 
and (3) by suitable differentiations and additions gives the 
following single equation in w for the equilibrium of cylindrical 
shells : 

By8w + it tjjjf +t v “(f It + 2T lv w 


Equation (4), which was first derived by Donnell (reference 11), 
was treated in reference 1. 


An alternative method for obtaining a single equation in w 
for the equilibrium of a cylindrical shell is to Bolve equation (3) 
for F and substitute the result into equation (2). This procedure 
can readily be carried out in the following manner. Differentiation 
of equation (3) twice with respect to x gives 




0 


( 5 ) 


The symbolic solution of equation (5) for 


d g F 

dx 2 


is 


d 2 F Jy4 d^w 


Substitution of this result into equation (2) gives 


2ft \?-4 , j. f zfiv ^ fjr ^ 

DV% + 3 S? + *(?* + 2T ar§; + spJ +£ ”° (6) 


Equation (6) is simply equation (4) modified by multiplication by 
the operator y^. m the present paper, equations (4) and (6) 
are referred to as Donnell* s equation and the modified equation, 
respectively. 


Advantages of Modified Equation 

One of the quickest and most convenient methods' for obtaining 
solutions of flat— plate buckling problems to any desired degree of 
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approximation uses a Fourier series type of expansion for the deflection 
surface w» Both Donnell's equation end the modified equation can 
"be solved "by this method in the case of "buckling problems involving 
curved plateshaving simply supported edges. 

As mentioned in the “Introduction," however, Donnell's equation 
is not well adapted to solution by Fourier series of problems 
involving the stability of shells with clamped edges . The cause of 
the trouble appears to be that the calculation of some of the high- 
order derivatives found in Donnell 's equation sometimes leads to 
divergent trigonometric series when the edges are clamped. The 
modified equation, however, is applicable to clamped-edge problems 
as well as to problems involving simply supported edges because 
lower-order derivatives are involved. 

In addition to its advantages in the solution of problems 
involving shells with clamped edge3, equation (6) has the additional 
advantage that each term, has a definite physical signif icance : 

The first term gives the restoring force per unit area of the deflected 
surface due to bending and twisting stiffnesses; the second term, 
gives the restoring force per unit area due to stretching stiffness; 
and the remaining terms give the deflecting forces per unit area 
due to applied loads. Because of these advantages, the modified 
equation was adopted for general use in references 2 to 7 • 

Both Donnell's equation and the modified equation result in 
the same critical stresses for simply supported cylindrical shells, 
and the two methods require essentially equivalent mathematical 
processes. (See appendix.) The characteristics of solutions by 
means of Donnell 'a equation in the case of simply supported shells 
(reference 1) - namely, the theoretical cylinder parameters, the 
simplicity of calculations and results, and the implied boundary 
conditions on u and v - are characteristics, also, of solutions 
by means of the modified equation. The same characteristics, except 
for a change in the implied boundary conditions on u and v, alBO 
apply to solutions of clamped-edge shell problems by means of the 
modified equation. This change is discussed in the section entitled 
'Boundary Conditions . " 


Solution of Modified Equation by Galerkin Method 

An approximate method of solving vibration and buckling problems 
closely paralleling that of Ritz was introduced in 1915 by Galerkin. 
(See, for example, references 12 and 13*) The main distinction between 
the Ritz and Galerkin methods is that the Ritz method begins with an 
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energy expression? whereas the Galerkin method, begins with an 
equation of equilibrium. The Galerkin method is readily adaptable 
to the solution of equation (6) and is now desoribed briefly. 

Let the equation of equilibrium be written 


Q(w) + p - 0 (T) 

where p is lateral pressure and. Q is scane operator in x and y 
which for the purposes of this paper is taken to be linear. According 
to the Galerkin method, the equation may be solved by expanding 
the unknown function w in terms of a suitable set of functions 
f i( x )sj(y)j each of which satisfies the boundary conditions but 
not in general the equation of equilibrium: 


w 



( x )gj(y) 


( 8 ) 


Substitution of this expression for w into equation (7) gives 
the following equation: 


/*~ a i jQ^iWejC?)] + p = o (9) 

i J ‘ 


Because the functions fi(x)gj(y) were chosen to satisfy the 

boundary conditions rather than the equation of equilibrium, 
equation (9) cannot, in general, be satisfied identically by any 
choice of the coefficients aij. These coefficients can be chosen, 

however, to assure the vanishing of certain weighted averages of 
the left-hand side of equation (9) . The weighting functions used 
in the Galerkin method are the original expansion functions, so 
that the following simultaneous equations for determining the 
coefficients aj. ^ are obtained : 



( 10 ) 


(m— 1 ,2,3, « • • > n=l ,2,3, • • • ) 


r 
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where 


B mni^ 


fm(x)sn(r)< Qj?i(x)gj(y)j + p > dx dy (11) 


The simultaneous set of linear algebraic equations in the unknown 
coefficients aij (equation (10)), obtained by using the original 

expansion functions as weighting functions, is ordinarily the same 
set which would be found by the Rltz method, if the same series 
expansion for w were used. A solution of any desired degree of 
accuracy may therefore be obtained by the Galerkin method. 


In applying the Galerkin method to equation (6) by use of 
Fourier series expansion for w, expressions of the type 



must be evaluated. The operator V the inverse of V 
simply Introduces into the denominator of each tepa of the series 
the expression that comes into the numerator if V ^ is applied. 
Thus, 



This result may readily be verified by applying the operator V 2 * 
to each side of equation (12) . 

In writing equation (12) the quantity V^f, as defined by 
the equation 
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vas tacitly assumed to be unique. The quantity actu all y is not , 
unique; any number of terms which vanish when operated upon by V 4 
could be added to the right— Land side of equation (12) . The 
omission of such terms makes the present analysis parallel to the 
analysis using Donnell’s equation (see reference 1) and implies 
certain boundary conditions on u and v, which are discussed 
in a subsequent section entitled "Boundary Conditions.", 


Deflection Functions 

Simply supported edges.- For simply supported cylindrical 
shells, the following series expansions for w may be used to 
represent the buckle deformation to any desired degree of accuracy 
(in these functions, x is the coordinate in the axial direction 
and y, the coordinate in the circumferential direction) : 

(1) Rectangular curved plate (axial dimension a and 
circumferential dimension b) 

00 CO 

V 

w = ^ y a^n sin sin (13) 

m~l n-1 


(2) Curved strip long in the axial direction (circumferential 
width b and axial wave length 2\) 

(a) Direct stresses only 


w = sin S V ani sin (l4) 

m-1 


(b) Shear stress with or without addition of direct 
stress 


00 0 , 

w * sin am sin + cos *> m sin (15) 

m= i m=l 


r 
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(3) Complete cylinder (length L and circumferential wave 
length 2X) 

(a) Direct stresses only 
00 

;Y 


w = sin sin 


m=l 


mjtx 

L 


(16) 


(h) Shear stress with or without addition of direct 
stress 




CO 


( 17 ) 


C lamped edges .— Probably the simplest method of treating 
cylindrical shells with clamped edges is to employ the expansions 
in equations (13) to ( 17 ) modified by substituting functions of the 
type 


<P m (x) = sin sin 


cos 


(m — l)nx 
a 


cos 


(m + l)itx 


( 18 ) 


wherever functions of the type sin ^22. appear, with a similar 

& 

substitution for functions of y (all terms involving summation 
subscripts m and n are thus changed; terms involving X, 

such as sin. remain unchanged). The functions cp^(x) form 

a complete set so that finite expansions for w of the type suggested 
for shells with clamped edges as well as those for shells with simply 
supported edges may be used to represent the buckle deformation 
to any desired degree of accuracy. 


Boundary Conditions 

Simply supported edges .— Appendix D of reference 1 shows that, 
if the buckling stress of a simply supported shell is found by means 
of the expansions for w given in the preceding section entitled 
"Deflection Function," the boundary or edge conditions implied for 
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the median-surface displacements u and v are zero displacement 
along each of the edges of a cylinder or curved panel and free 
displacement normal to each edge . (Although the proof given used 
equation (4), the proof could equally well have "been "based on 
equation (6),) 


The "boundary conditions for simple support may thus be -written, 
at a curved edge (x = constant) , 


d 2 * 5% 

Sx 2 = T = Sy2 


0 


( 19 ) 


and, at a straight edge 


(y = constant). 


S^w 5 2 F 


( 20 ) 


Clamped edges . - By a method similar to that in appendix D 
of reference 1, solutions using the functions suggested in the 
preceding section for the treatment of clamped edges can be shown 
to correspond to the boundary conditions - zero displacement normal 
to an edge and free displacement along an edge . 

The boundary conditions for clamped edges therefore become, 
at a curved edge (x = constant), 


w 'H” u = a^-° 


( 21 ) 


and, at a straight edge (y = constant). 


3w __ d^F n 
dy dy 2 


( 22 ) 


Discussion . - As mentioned in reference 1, the boundary conditions 
implied for u and v in the case of simply supported edges are 
appropriate for cylinders or panels bounded by light bu lithe ads or 
deep stiffeners, which are stiff in their own planes but may be 
readily warped. out of their planes. 
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The boundary -conditions on u and v appropriate for a 
clamped edge would seem to be zero displacement normal to the edge 
and zero, rather than free, displacement along the edge. Comparison 
of critical stresses for shells with clamped edges found by the 
method in the present paper with critical stresses found by the 
method in references 9 and lij-, giving boundary conditions u » v » 0, 
however, indicates that the imposition of the added requirement 
of zero displacement along, the edge ordinarily has very little 
effect on the critical stresses. 

A less satisfactory method for solving problems concerning 
shells with clanged edges involves the use of functions of the type 


1 

m 


sin fflUS _ _JL_ 
sxn a m + ;2 


sin 


(m + 2) its 
a 


instead of those described by equation (l8) . In this method, the 
functions used are those for single support taken in such combinations 
that the edge slope is zero. Use of such functions leads to the 
same boundary conditions on u and v as were described for Bimply 
supported edges; at the edge y = constant, for example, the boundary 
conditions become 


w 


& 

3y 


u 


a g F 

5x 2 


o 


( 23 ) 


The use of these functions to represent shells with clamped edges 
is not recommended, however, for the following reasons: The associated 
boundary conditions seem to be artificial and unlikely to be reproduced 
even approximately in actual construction.; the method leads in some 
oases to solutions that differ considerably from the solution for 
ideal clamped-edge conditions in which u = v ■= 0; and the solutions 
obtained generally converge rather poorly. 


APPLICATIONS AND DISCUSSION 


Among the more difficult shell -stability problems to treat 
theoretically are those which Involve shear stresses. In fact, 
until 1934, when Donnell’s paper on critical shear stress of a 
cylinder in torsion was published (reference 11), such problems 
wears generally regarded as impracticable to solve . In order to 
Illustrate the type of solution to be found by the method of analysis 
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just outlined and the effect of boundary conditions on critical 
stresses, the results obtained for a number of shell stability 
problems involving shear stresses are reproduced and discussed 
briefly here. The problems treated are summarized in table I. 

C ritical shear stress of long curved strip. - The critical 
shear stress for a long plate with transverse curvature is given 
by the equation 


T 


cr 



where 3% is a dimensionless coefficient, the value of which 
depends upon the dimensions of the strip, Poisson's ratio for 
the material, and the type of edge support. In figure 1 (fig. 1 
of reference 2) the shear -stress coefficient k s is given for 
plates with simply supported edges and with clamped edges. This 
solution for simply supported edges coincides with that given by 
Kromm (reference 15). 

As indicated in the previous section entitled 'Poundary 
Conditions, " the solution corresponding to the boundary conditions 
of equation (23) (dashed curve of fig. 1) is poorly convergent and 
deviates appreciably from the results for completely fixed edges. 

Figure 1 shows this poor convergence in the limiting case of a flat 
plate, for which the critical stress is independent of boundary 
conditions on u and v. Even a tenth -order determinant led to 
a result that is 7 percent above the true solution-, whereas the 
result using a fourth-order determinant obtained with the deflection 
functions recommended for clamped edgeB is only 1 percent above. 

In figure 2 (fig. 2 of reference 2} the solutions given in 
figure 1 are compared with the results given by Leggett (reference 9) 
for simply supported and clamped edges with u = v = 0 at each 
edge. Throughout the rang© for which they are given, Leggett's 
results for clamped edges differ only slightly from those of the 
present paper. On the other hand, the previously mentioned 
discrepancy between the results for completely fixed edges 
(u = v = 0) and those for the boundary conditions of equation (23) 
(dashed curve) may be inferred from this figure to be considerable 
for large values of Z. A minimum measure of this discrepancy is 
the distance between the clamped -edge curves for v = 0 and for 
u = 0 in figure 2, since Leggett’s curve must always lie above the 
curve for v = 0. 

The reason for the marked increase in buckling stress of simply 
supported curved strips when the edges are restrained against 
circumferential displacement during buckling is discussed In reference 2. 
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C ritical shear stress of cylinder in torsion . - The critical, 
shear stress of a cylinder subjected to torsion is given by the 
equation 


i£l> 


In figure 3 (fig* 1 of reference 3) the values of k s are given 
for cylinders with simply supported edges (boundary conditions 
of equation ( 19 )) and cylinders with clamped edges (boundary 
conditions of equation (21)). At high values of Z, the values 
of k s for thick cylinders are given by special curves for various 

values of ^ \J 1 - |j2, as discussed in reference 1. At values 

of Z greater than about 100 only a small increase in buckling 
stress is caused by clamping the edges. The results indicated 
in figure 3 are in very close agreement with Donnell's results for 
the same problem, except in the range 5 < Z < 500 wheye the 
somewhat lower curves of the present paper represent a more accurate 
solution. 

Reference 1 shows that boundary conditions imposed upon u 
and v at the curved edges of a panel or cylinder have an almost 
insignificant effect on the buckling stresses, whereas conditions 
imposed on v at the straight edges may be quite important. 
Comparison of figure 1, in which boundary conditions on straight 
edges are considered, with figure 3 , in which conditions on curved 
edges are considered. Indicates that a similar situation exists 
with respect to restraint against edge rotation. 

Critical shear stress of curved panel . - The values of k s 
giving the critical shear stresses of simply supported curved 
rectangular panels are given in figures 4 and 5 (figs. 1 and 2, 
respectively, of reference 4) . The corresponding boundary conditions 
on u and v are aero displacement parallel to the edges and free 
warping normal to the edges . Figure 4 indicates that as the 
curvature parameter Z increases, the critical shear stresses of 
panels having a circumferential dimension greater than the 
axial dimension approach those for a cylinder. Figure 5 indicates 
that, as the curvature parameter Z increases, the critical shear 
stresses for panels having an axial dimension greater than the 
circumferential dimension deviate more and more from the critical 
shear stress for an infinitely long curved plate. Reference 4 
shows that the reason for this deviation in figure 5 is that at 
high curvatures the buckling stresses of -these panels, as well as 
those of figure 4, approach those of the cylinder obtained by 
extending the circumferential dimensions of the panels. 
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The effects of "boundary conditions in the limiting cases of 
infinitely long curved strips (fig,. 1) and of complete cylinders 
(fig. 3) suggest that the curves of figure 4 sire substantially 
independent of edge restraint at large values of Z "but that the 
curves of figure 5 would "be considerably affected by a change in 
edge restraint. 

Long curved strips tinder combined shear and direct axial stress . 
Reference $ shows that the theoretical interaction curve fbr a 
long curved strip under combined shear stress and direct axial stress 
is approximately parabolic when the edges are either simply supported 
or clamped, regardless of the value of. Z. This parabola is given 
by the formula 


R 0 2 + R x = 1 


where R a and Rx are the shear -stress and compressive -stress 
ratios, respectively. 

At high values of Z curved strips, like cylinders, buckle 
at compressive stresses considerably below the theoretical critical 
stresses. In order to take this condition into account, certain 
modifications in the theoretical results are proposed in reference 5 
for use in design. 

Cylinders under combined shear and direct axial stress . - 
The theoretically determined combinations of shear stress and direct 
axial stress which cause a cylinder with simply supported end 
clamped edges to buckle are shewn in figure 6 (fig. 1 of reference 6) . 
Considerable variation in the shape of the interaction curves 
occurs for low values of Z. For high values of Z the interaction 
curves for either simply supported or clamped edges are similar to 
the curve for Z = 30- 

Because cylinders actually buckle at a small fraction of their 
theoretical critical compressive stress, the theoretical interaction 
curves of figure 6 cannot be expected to be in satisfactory agreement 
with experiment when a very appreciable amount of compression is 
present. For semi empirical curves and a check of available tost 
data, see reference 6. 
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CONCLUDING- KEMAKKB 


A previous investigation showed how Fourier series type 
solutions of Donnell’s equation can "be used to simplify greatly 
the stability analysis of thin cylindrical shells with simply 
supported edges . The present paper shows that the restriction to 
simply supported edges can be removed by the introduction of a 
new equation which is equivalent to Donnell 's equation but is better 
adapted to solution by Fourier series. This modified equation can 
be solved for the buckling stresses of curved sheet having either 
simply supported or clamped edges by established methods essentially 
equivalent to those in use for flat sheet. This approach permits 
a simple and straightforward solution to be given for a number of 
problems previously considered rather formidable. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., March 20, 19^7 
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APPENDIX 

COMPARISON OF RESULTS OBTAINED BY USING DONNELL'S EQUATION 
AND TEE MODIFIED EQUATION IN THE STABILITY ANALYSIS 
OF SIMPLY SUPPORTED CURVED PANELS 
Solution of Donnell's Equation 


Donnell's equation expressing the equilibrium of a curved 
panel under median-surface stresses can be written in general form 
as • 


8 Et 

DVV+ i2£5 + 


o x tV 


kcfv 
3x 2 


+ 

ox dy J 


sA 5% „ 0 

5y 2 


(Al) 


where x is the axial coordinate and y the circumferential 
coordinate. Division of equation (Al) by D and the introduction 
of the dimensionless stress coefficients k x , k-j, and k s , and 
the curvature parameter Z results in the following equation: 


,_8 ^ 12 Z 2 

v W + — r— — r 

b^ cte 5 


1 r-~r 


k 8“ W 




a ^ 2 + 21t V v ax a y 






8y 2 


0 (A2) 


where 


and 


b 2 t 


** - %“55 


V b 2 t 

** 7 Zd 


k - a — 

y ' Vi 


z = R /i - H £ 


r 
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Equation (A2) can be represented by 


Ql(w) =* 0 


(A3) 


where Q^_ is defined as the operator 




The equation of equilibrium (equation (A3) ) is solved by using 
the Galerkin mettled as described in the section entitled "Theory. " 

In applying this method the unknown deflection w 1 b represented in 
terms of a set of functions (see equation (8)), each of which 
satisfies the boundary conditions but not in general the equation of 
equilibrium. A suitable set of functions of thlB type, which 
satisfies the boundary conditions for simple support. Is 


CO CO- 

% ain 5S sin (A4) 

m=l n=l 


where the origin is taken at a corner of the plate Substituting 
in equations (10) and (ll) 



fm( x ) 


*= sin 


BUtX 


eh(y) = sin ^ 


and 


Q = Qq 


and performing the integration over the whole plate (limits x = 0, a 
y = 0, b) gives the set of equations 
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®Hn 


(A5) 


where m=l,2,3, «.*, n»l,2,3> • and P and q take only those 
values for which m ± .p and n ± p are odd numbers . 

Equation (A5) represents an infinite set of homogeneous linear 
equations involving the -unknown deflection coefficients a^ ^ . In 

order for the deflection coefficients to have values other than 
zero, that is, in order for the panel to buckle, the determinant of 
the coefficients of the unknown deflection coefficients a., . 

must vanish. This determinant can be factored into two subdeterminmts, 
one involving the unknown deflection coefficients a^ j for which 

i * j la odd and the other involving those coefficients for which 
i i J ie even. Buckling occurs, therefore, when either of the two 
subdetflirminants vanishes. Only the buckling criterion Involving 
the even subdeterminant is treated here . This criterion is 


- k x m2 (^2 + n 2 - k^n 2 §§ (f + n 2 20 


32kr s a3 

it 2 b^ 


P=1 1= 


(p 2 + q 2 mnpq 

(5T— 2K -jrr^) - 0 


£ + $) k * 


r 
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a ll a 13 a 22 

®3l 

®33 

m=l, n=l 

bln 0 +£ (4-+4^j 

2 

0 

0 

m=l, n=3 
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where 



Division of each column of the 
by the proper 


determinant in equation (A 6 ) 



a 0 (A 6) 
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gives the simplified equation 
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where 
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The vanishing of this determinant is the criterion for the symmetrical 
bucJcling of the shell. The same buckling criterion results from 
the use of the modified equation, as is shown in the following 
section. 


Solution of Modified Equation 

The modified equation expressing the equilibrium of a curved 
panel under median-surface stresses in general form is 
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(A8) 


Division of e a nation (A8) tv D and simplification of the result 
gives the following equation: 


V^w 


w + t, «g + 2k„ i|25k_ + te, £ . 0 (A9) 
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Equation (A9) can he represented ty 


QgCw) = 0 


(A10) 


where Qg is defined as the operator 


V* 


+ 




« 2 82 




* 2 82 , , £ 82 _ 
b 2 8x 8y ^ t 2 Sy 2 


By use of the Galerkin method and ty use of the expression for v 
given in equation (A4 ) , the following set of equations analogous 
to equations (A5) are obtained 


i mn 
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(m 2 - p 2 )(n 2 - q 2 ) 

P=1 1=1 (All) 


where m=l,2,3, n=l,2,3, •••, and p and q take only those 

values Buch that mt j and n ± q are odd numbers. 
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As In the case of the solution of Donnell's equation, the 
stability determinant representing equations (All) can he factored 
into an even and an odd suhdetsrminant. The even one is 
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m=3. 

n=3 

0 

0 

+.3S 

0 

Noo 


2*5 


=> 0 (A12) 


The stability determinant (equation (A12)) obtained from the 
modified equation is identical with the simplified stability 
determinant (equation (A7)) obtained by use of Donnell's equation. 
This identity holds for the odd as well as the even determinants. 


Although the stability determinants obtained by use of the two 
equations are identical and yield identical buckling loads, the 
determinant in equation (A7) consists of the coefficients of 



whereas the determinant in equation (A12) 


consists of the coefficients of Accordingly, although the 

buckling loads found by the two methods are the same, the buckle 
patterns are different. Of the two buckle patterns the one found 
by the use of the modified equation is believed to be correct. 
This conclusion has been verified for the limiting case of a 
flat plate (Z=0) . 
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Figure 2.- Comparison of Leggett’s solutions with present solutions for 
critieal-shear-stress coefficients of a long curved strip. (Fig. 2 of 
reference 2.) 
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Figure 5.- Critical-shear-stress coefficients of simply supported 
curved panels having axial dimension greater than circumferential 
dimension. (Dashed curve 'estimated.) (Fig. 2 of reference 4.) 
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(b) Clamped edges. 


Figure 6.- 
stress 


Critical combinations of shear-stress and direct-axial 
coefficients for cylinders. (Fig. 1 of reference 6.) 














